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Abstract 

The asymptotic form of the energy density for a gas of particles surrounding 
a sphere of mass M and radius R is studied using Einstein's equations. It 
is shown that if the pressure of the gas p varies linearly with the energy 
density p for small p, then p ~ 1/r 2 for large r. 
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In the Schwarzchild solution of Einstein's equations for a static, spheri- 
cal geometry [1] the spacetime surrounding the sphere of mass M and radius 
R is taken to be empty and free of particles. The real world is certainly not 
so clean and simple, however and a more realistic model would also include 
a gas of particles which surround the sphere. With the addition of these 
particles the question then becomes what affect, if any, they will have on 
the geometry of the spacetime. At first glance this would seem to be an 
impossible question to answer. Not only must an equation of state for the 
particles be known, this equation of state must also have been calculated in 
the presence of the very gravitational field that it itself determines. What 
we shall find, however, is that due to the restrictions that Einstein's equa- 
tions themselves put on the energy density of the particles, the asymptotic 
behavior of this spacetime will be surprisingly simple to determine. 

In this paper we shall study some of the asymptotic properties of a 
static, spherical geometry filled with a gas of particles. To be specific, the 
system we shall be considering consists of a sphere with mass M and radius 
R which is surrounded by an arbitrary gas of particles. We shall not need 
to specify the type of particles surrounding the sphere, nor shall we need 
to fix their temperature or density. All that we shall require is that they 
have an equation of state which has certain properties. For convenience we 
have also assumed that the sphere has not undergone complete gravitational 
collapse into a black hole and that the system as a whole is in thermodynamic 
equilibrium. Moreover, unlike the usual treatment of matter in general 
relativity, we shall not a priori confine the particles to be within a sphere 
of any fixed radius, but will instead let the system itself determine how far 
into the spacetime the particles will spread. 

Although we have used the word "gas" to describe the particles, this is 
only a matter of convenience. The particles are not only allowed to interact 
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with gravity, thereby determining the geometry of the spacetime, but also 
with each other. Since the most stringent constraint on the energy density, 
and thus the geometry of the spacetime, comes not from the equation of 
state of the particles, but rather from Einstein's equations themselves, all 
that we need require is that the equation of state satisfy a few physically 
reasonable requirements. In the end we shall find that if the particles in- 
teract among themselves to any significant extent, then the spacetime will 
be asymptotically flat. The energy density of the particles will die off faster 
than 1/r 2 and we may characterize the particles as being confined within a 
sphere of a certain radius which will require a specific choice of the equation 
of state to determine precisely. If the particles are very weakly interacting, 
on the other hand, then we shall find that the spacetime is not asymptoti- 
cally flat. The energy density will die off as 1/r 2 for large r and we cannot 
characterize the particles as being confined within a sphere of any fixed ra- 
dius. The affect of the particles on the geometry of the spacetime is instead 
all pervasive. 

We begin with the most general form of the metric for a static, spherical 
geometry [2] 

ds 2 = -fdt 2 + hdr 2 + r 2 d6 2 + r 2 sin 2 Odcf 2 , 

where / and h are unknown functions of r only. We are only interested in 
the structure of the spacetime for r > R and shall always assume that this 
holds. As usual, we shall take the average energy momentum tensor for the 
particles to be of the form (T^) = pu^u^ +p{g lJLV +u^u v ) where n M is a unit 
velocity vector which lies in the direction of the timelike Killing vector for 
the system, while p and p are the energy density and pressure, respectively, 
of the particles surrounding the sphere. The two Einstein's equations we 
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shall find of use are [2] 



87vp -JL + 1 ( i _ l\ 
rh 2 r 2 \ h J ' 

rfh r 2 \ h J 



(la) 



(lb) 



where the primes denote derivatives with respect to r and we are using units 
in which G = c = 1. From the conservation equation V At (T At!y ) = we also 
have the Tolman-Oppenheimer-Volkoff (TOV) equation [2] for hydrostatic 
equilibrium 



At this point we shall assume that we are given an equation of state for 
the particles which we shall write as 



where w(p) is a functional of p. It is further assumed that this equation 
of state includes the rest mass of the particles and was calculated in the 
presence of the gravitational field. As such both p and p must also be 
functions of r. We next make the anzatz that w(p) itself has no explicit 
r dependence, but is instead dependent on r only through its dependence 
on p. We justify this anzatz with the following. First, we know that in 
Minkowski spacetime the equation of state may always be written in the 
form given in eq. (3). As the system is in thermodynamic equilibrium, and 
as there are no external fields present, any r dependence in w(p) must then 
be due to the curvature of the spacetime. Second, log / may, under certain 
circumstances, be interpreted as being twice the gravitation potential and 
is therefore a measure of the "gravitational force" on the particles. Third, 
we note that all explicit r dependence in the TOV equation can be canceled 
and one can instead consider it as a differential equation determining p in 




(2) 



p = w(p)p, 



(3) 
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terms of /. In fact, any dependence of p on r is due to /. Finally, the 
geometry of the spacetime, namely / and h, is ultimately determined by 
the energy density and pressure of the particles. Since the equation of state 
on curved spacetimes must also reduce to the equation of state calculated 
in Minkowski spacetimes in the limit where / — > 1 and h — > 1, we would 
therefore expect w(p) to depend on r only implicitly through p. 

The equation of state will in general be very complicated, if it can be 
calculated at all. Fortunately, we shall not need any specific form of eq. (3), 
but rather that it have a couple of physically reasonable properties which 
we would expect from any physical system. First, the pressure p must be 
finite as p — > 0, meaning that the gas of particles cannot become infinitely 
stiff as the energy density of the particles is reduced to zero. Actually, we 
shall need the more stringent requirement that 

lim w(p) = Wq , 

exists and is a finite. This means, in particular, that for small p, p ~ p 1+n for 
n > 0. We justify this requirement with the observation that if the particles 
were consisted only of photons, then the pressure must be one third the 
energy density. This holds even in curved spacetimes. Since the photons 
do not interact among themselves, we have at least one example of a gas 
of non-interacting particles for which the pressure varies linearly with the 
energy density. We would therefore expect that if the particles do interact 
with one another, then the pressure should vary at least linearly with the 
energy density, if not stronger. We would not expect the pressure to have a 
weaker power law dependence on the energy density than a linear one. 

With the Einstein's equations, and the TOV equation, we have now 
three differential equations determining p, /, and h. (Since the equation of 
state is given, p is determined as soon as p is.) We shall then, presumeably, 
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need a set of three initial conditions to determine the system completely. 
The initial condition for p is straightforward; we need only take take its value 
at the surface of the sphere pn = p(R). The initial conditions for / and h, 
on the other hand, are much more difficult to determine. Fortunately, we 
are only interested in the asymptotic nature of the spacetime and we shall 
later find that the asymptotic forms of both p and h are independent of 
their initial conditions. Whatever initial condition we choose for h is thus 
immaterial for our purposes. /, on the other hand, does depend on initial 
conditions; not only on its own, but also on pr as well. Since, however, we 
may always rescale the time coordinate, its dependence on its own initial 
condition is not extremely relevant and we may formally take it to be f(R). 

Because the equation of state may be quite complicated, the solution 
of eqs. (1) for all r is not trivial to find. In, however, the asymptotic limit 
where r — > oo the situation simplifies dramatically. To see how this happens, 
we integrate eq. (la) to formally give 



for r > R where mo is an integration constant. Since p is not known 
explicitly, this would seem to be of little use. Note, however, that h > 
for all r > R. Consequently, p — > as r — > oo. In fact, p must die off at 
least as fast as 1/r 2 . Thus, because we are only interested in the asymptotic 
behavior of the spacetime, only the behavior of the equation of state when 
p — > 0, namely wo, will be relevant in our analysis. Keeping this in mind, 
we next make use of the equation of state in eq. (2) to obtain 






We then define 



A 



P = 
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and note that A either vanishes as r — > oo or else approaches a constant 
value. The important point is that A must be finite as r — > oo. Then 
defining h~ x = 1 — 2K, eqs. (1) become 

' iJK +K, (5a) 



dy 

w {l-2K) ( IdA \ r , 

Aw = ^ -r— 2 ) — K , 56 

l + w \A dy J 

for large y where y = log (r/r ) for some tq. Solutions to these equations 
have two different behaviors depending upon whether or not w$ vanishes. 
Case 1: wq = 

From eq. (56), K — > as r — > oo. Consequently from eq. (5a), A — > 
as well. This means that the energy density, and consequently the pressure, 
must decrease faster than 1/r 2 for large r. Then from eq. (4), / — > Ck 
where Ck is a constant which can be set to 1 by suitably rescaling the time 
coordinate. Consequently, if wq = 0, then / — > 1 and a — > 1 as r — > oo and 
the spacetime is asymptotically flat. 
Case 2: w ^ 

This case is much more interesting. We write eqs. (5) as a set of two 
coupled, non-linear differential equations 

^=A-K, (6a) 
dy 

1+W " A ^ roo+ fi±^W-^j. (66) 



a 1 ?/ wo 1- 2K \ \ 1 + wo ) 1 + 

Fortunately, these two equations are autonomous, meaning they have no 
explicit y dependence. Obtaining asymptotic solutions to eqs. (6) is then 
straightforward and involves looking for fixed points (A a ,K a ) of eqs. (6) 
where the derivatives of both A and K vanish. (See [3]. We caution the 
reader that, depending on the literature, the term "critical point" is often 
used instead of "fixed point" .) From eqs. (6) one of these fixed points occurs 
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at A a = = K a , which is the wq = case once again. A second fixed point 
occurs at 

A -K - 2W ° 



(1 + w ) 2 + 4w ' 
If we next expand eqs. (6) about this fixed point, then 

where 5 A = A — A a and 5K = K — K a . The asymptotic behavior of the 
solutions to eqs. (6) depends on the eigenvalues A± of this matrix. Writting 
A± = — r\ ± i<p, we find that 

1/ 2w \ l(7 + 42w Q -w 2 ) 1 / 2 

V = o 1 + T~, ' t P~ 



2 \ 1 + wqJ 2 l + w 

Since p > for all r, A a must be positive. This condition, combined with 
the requirement that h~ x = 1 — 2K a > for all r, means that wq > 0. 
Moreover, because p < p/3, wq < 1/3, and it is then straightforward to see 
that ip will always be a real number. Consequently, this fixed point is stable 
and all solutions to eqs. (6) will eventually spiral counterclockwise into this 
fixed point in the A-K plane no matter what their initial conditions were 
originally. We can see this explicitly by solving eq. (7) to give 

• I 1 r 

sin ip log — 

r , 

T 

if COS I if log — 

r 





+ (1 - rj) sin \ip\og — 
V r 

for large r. A, and r$ are constants which depend on the specific equation 
of state and initial condition for p and p. Notice also that because < 
w < 1/3, < A a < 3/14, and 1/2 < 77 < 3/4 while V7 '/2 < </? < 
Solutions to eqs. (6) therefore approaches this fixed point very slowly; the 
fastest being r~ 3 / 4 . 



From eq. (4) we can solve for / in terms of p. In the asymptotic limit 
we find that / rs kr q where q = 4w /(l + wq), 



and wr = w(pft). Since < w < 1/3, < q < 1. Consequently, for large r 
the metric becomes 



and we can see explicitly that this spacetime is not asymptotically flat. 

To summerize, we have shown that if the particles interact among them- 
selves to such an extent that w = 0, then a static, spherical spacetime con- 
taining these particles is asymptotically flat, as expected. Moreover, we see 
that the energy density of these particles must necessarily decrease faster 
than 1/r 2 ; most probably very much faster. As the Schwarzchild geometry 
is also asymptotically flat and involves a mass density which is confined 
within a sphere of definite radius, we may therefore characterize these par- 
ticles as effectively being confined within a sphere of a certain radius which 
will require a specific equation of state to determine. The affect of these 
particles on the curvature of the spacetime, like the affect of the mass M in 
the Schwarzchild solution, eventually dies away and the spacetime asymp- 
totically approaches the Minkowski spacetime. When, on the other hand, 
the particles interact with each other so weakly that w ^ 0, then their 
energy density decreases as 1/r 2 and their affect on the curvature of the 
spacetime is correspondingly long range. They cannot be characterized as 
being contained within a sphere of any definite radius and are instead spread 
throughout the spacetime. In fact, we find that / ~ r q while h approaches 




ds 2 = -kr q dt 2 + 



1 



dr 2 + r 2 d6 2 + r 2 sin 2 Odcf) 2 , 



1-2K, 



a 
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a constant as r — > oo and a spacetime filled with these weakly interacting 
particles is not asymptotically flat. 

What is of even greater interest is the universal nature of the energy 
density in the asymptotic limit 

A a c 4 

Pa = 4^G ' (8) 

where we have replaced the correct factors of c and G. (The subscript 
a denotes the asymptotic limit.) Notice that the form of p a is the same 
irrespective of the mass and radius of the sphere, irrespective of any "initial 
condition" for p at r = R, and irrespective of any form that the equation of 
state may take so long as w ^ 0. In fact, the only property of the particles 
that it does depend upon is wq, a dimensionless number. From dimensional 
arguments, if the particles surrounding the sphere are massless, then wo 
must simply be a number. It cannot even depend on the temperature of the 
gas, since, aside from the Planck mass, there is no other length scale one 
can use to construct a wq. If, on the other hand, the particles have a mass 
m, then wo is either once again a number, or else is a function of the ratio 
T/m, where T is the temperature of the gas. 

Notice that although this universality extends to h, it does not extend 
to include /. For large r, / ~ kr q and although q depends only on wo, 
k depends not only on f(R) but also on p^. Consequently, the motion of 
test particles in this spacetime will always be dependent the specific choice 
of initial conditions for / and p, and thus on the detail properties of the 
particles and of the spherical mass. 

The question remains as to how one may go about calculating wq- 
For certain cases this is trivial. Suppose T^ v is the energy momentum 
operator for the particles such that its equilibrium average is (T Ml/ ) = 
pu^Uy + p(g^ u + u^u v ). If this energy momentum operator is traceless, 
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as it is for a gas of pure photons, then (T^) = and p = p/3. Conse- 
quently, any system which has a traceless energy momentum operator has 
a w(p) = wq = 1/3. Calculating wo for other systems, on the other hand, 
is a much more formidable task, although at first glance not an impossible 
one. Because wq is determined when the energy density of the particles is 
vanishingly small, this suggests that when calculating wq one may, as a first 
approximation, neglect the affect of the particles on the curvature of the 
spacetime and may instead treat them as test particles. One can then in 
principle use the vacuum solutions of Einstein's equations as a background 
field and calculate the equation of state for the system using methods de- 
scribed in [4] . 

We end this paper with a couple of observations. First, the 1/r 2 be- 
havior in eq. (8) is precisely what one would expect for the energy density 
of "dark matter" based on Newtonian gravity [5] — [7] . Unfortunately, when 
one has gone to such a large r such that eq. (8) holds, one also finds that 
the spacetime no longer close to being Minkowskian. It is instead extremely 
curved (see [8]). Second, the form of p a is quite peculiar in that it contains 
no explicit length scale for r. The only length scale which can be constructed 
solely from universal constants, however, is the Planck length = hG/c 3 . 
If we use it, then we can write eq. (8) as 



where p p i = c 7 /(hG 2 ) is the Planck energy density. Although p p i is very 
large, eq. (9) is valid only at very large r. Since l p i ~ 10 _33 cm, this ensures 
that the size of p will always be physically reasonable. 
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